Magnetism in 2D BNi_a;Oa; and Bi_a;Sia;N: 
polarized itinerant and local electrons 
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We use density functional theory based first-principles methods to study the magnetism in a 2D 
hexagonal BN sheet induced by the different concentrations of oxygen and silicon atoms substituting 
for nitrogen (On) and boron (SIb) respectively. We demonstrate the possible formation of three 
distinct phases based on the magnetization energy calculated self-consistently for the ferromagnetic 
(MEfm) and antiferromagnetic (MEafm) states, i.e. the paramagnetic phase with MEfm=MEafm, 
the ferromagnetic phase with MEpm>MEafm and finally the polarized itinerant electrons with 
finite MEpM but zero MEafm- While the On system was found to exist in all three phases, no 
tendency towards the formation of the polarized itinerant electrons was observed for the Sis system 
though the existence of the other two phases was ascertained. The different behavior of these two 
systems is associated with the diverse features in the magnetization energy as a function of the 
oxygen and silicon concentrations. Finally, the robustness of the polarized itinerant electron phase 
is also discussed with respect to the O substitute atom distributions and the applied strains to the 
system. 

PACS numbers: 75.75.+a, 71.10.Pm, 71.10.Ca, 61.72.-y 
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I. INTRODUCTION 

The theoretical model of the two-dimensional (2D) 
homogeneous electron gas (HEG), which considers the 
electrons moving in a 2D uniform positive charge back- 
ground, serves as an important model system for study- 
ing the fundamental many-electron behavior in 2D sys- 
tems. These physical properties are closely connected 
to the operations of conduction electrons in the layered 
semiconductor devices as well as the recently synthesized 
2D periodic systems [ll, The most creditable method 
considered presently for studying the 2D HEG is the 
quantum Monte Carlo (QMC) method. The recent 
QMC studies found no region of stability for a ferromag- 
netic fluid in the 2D HEG 6]. The phases were found to 
transit directly from a polarized Wigner crystalQ to the 
normal fluid (paramagnetic state). However, it would be 
interesting to explore the parallel phase diagram for real 
materials, i.e. systems of electrons moving in a neutral- 
izing background formed by a 2D lattice of nuclei, and 
study the corresponding magnetic property as a function 
of the electron density. 

Inspired by the exciting experimental works on the 
discovery and synthesization of hexagonal boron nitride 
(h-BN) in low-dimensional structures [l|, 0, H, as well 
as the possible formation of magnetism in these systems 
through creations of defectspij jlll . here we present Den- 
sity Functional Theory (DFT)[12ll based first-principles 
studies on the magnetic properties of a 2D hexagonal 
BN sheet induced by the oxygen and silicon atoms sub- 
stituting for the N and B atoms respectively (denoted as 
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On and SIb hereafter), i.e. BNi_a;02: and Bi-ajSi^jN. We 
demonstrate the possible formation of polarized itiner- 
ant electrons in the On system. The stable phases are 
identified as, in increasing substitute atom concentration 
X, the paramagnetic phase, the ferromagnetic phase and 
finally the polarized itinerant electrons. In addition, the 
SIb system which though supports the formation of ferro- 
magnetic phase displays no tendency to the stabilization 
of the polarized itinerant electrons, i.e. resembling the 
QMC result for a 2D HEG. 



II. COMPUTATION METHOD 

The DFT calculations employed in the present st udy 
use the generalized gradient approximation (GGA)|13l| 
for the exchange-correlation energy functional and the 
projector augmented- wave (PAW) method[l3,[ll| for de- 
scribing the interaction between core and valence elec- 
trons as implemented in the VASP^] code. To sim- 
ulate the systems with the substitute atom concentra- 
tion X ranging from 1/55 to 1/4, the supercells with dis- 
tance (denoted as dnn) of ~17.5A, 13. SA, lOA, 7.5A and 
5A between nearest-neighbor substitute atoms were con- 
structed. A vacuum distance larger than 12A was used 
in the calculations to remove the interaction between lay- 
ers. The numerical convergence was accomplished by us- 
ing k-points meshes of (7 7 1) generated from Monkhost- 
Pack method [l3| for the supercell consisting of 2 x 2 unit 
cells (and the similar density of k-points meshes for the 
supercells of other sizes), and the kinetic-energy cut-off 
of 400eV (with tests at 500 eV) for expanding the sin- 
gle electron Kohn-Sham wavefunctions in plane- wave ba- 
sis. The atomic forces calculated by Hellmann-Feynman 
Tr^v.^th-eoremfisjl were relaxed to less than 0.02eV/A in all 
calculations. Furthermore, the volume and shape were 
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FIG. 1: The magnetization energy MEfm (solid lines) and 
MEafm (dashed lines) per substitute atom with respect to 
dnn as well as the substitute atom concentration x. The circles 
and triangles are for the systems of On and SIb respectively. 

also allowed to relax for the systems constructed by 2 x 2 
(x = 1/4 and dnn^SA) and 3x3 (a; = 1/9 and dnn^T.SA) 
supercells. We shall emphasize that the relaxed struc- 
tures for On and SIb are found to likely preserve its 
three-fold symmetry in the hexagonal lattice at all stud- 
ied substitution concentrations |19|. When switching on 
the spin-polarized calculations, the atomic forces were 
further relaxed. 



III. RESULTS 

The calculated magnetizations for both On and SIb 
stay as IfiB per substitute atom throughout the systems 
of fractional x's studied in this work[2(),]. The magnetic 
moment found in both systems can be attributed to the 
donated electrons accommodating at the bottom of the 
conduction band due to the presence of the substitute 
atoms which is confirmed by the corresponding density 
of states (DOS) of these systems. The magnetization en- 
ergies (ME) for the ferromagnetic and antiferromagnetic 
configurations are defined as the energy gains of these 
configurations (EpM and Eafm) with respect to those ob- 
tained from the non-spin-polarized calculations (Ensp), 
i.e. MEfm=Ensp-Efm and MEafm =Ensp -Eafm re- 
spectively. (All the magnetic related quantities discussed 
in this work are given in the unit of per substitute atom.) 
Assuming that the interaction among the local mag- 
netic moments in these magnetic systems can be approx- 
imately described by the nearest-neighbor Heisenberg 
Hamiltonian, H = — X)i>j JSiSj , then the interaction 
strength as well as the magnetic coupling types, i.e. fer- 
romagnetism and anti-ferromagnetism, can be extracted 
from the exchange energy J = ^(MEfm— MEafm)- Note 



that the initial magnetic configurations in calculating 
Eafm were set up with antiparallel local moments of 
equal magnitude and then followed by the spin-density 
relaxation in the self-consistent calculations. In FIGH] 
we present the calculated magnetization energy of On 
and SIb as a function of the substitute atom concen- 
tration. At low concentrations of x < 1/55 both sys- 
tems are stabilized at paramagnetic phase, i.e. J — 
but with finite MEfm=MEafm of around 50meV for On 
and 140meV for SIb. The presence of substitute atoms in 
these concentrations therefore introduces local moments 
into these systems yet no interactions among them, i.e. 
isolated local moments. These moments were identified 
as mostly Pz-^ike and their locality was manifested from 
the narrow band-width bands splitting from the original 
conduction band of the BN sheet. The spin density for 
SIb distributes mostly around the Si atom sites, while 
that for On is considerably more extensive, i.e. up to the 
third nearest neighboring B atoms to the O atom. 

As the substitute atom concentration is increased, both 
systems develop into the ferromagnetic phase with finite 
J. The values of 4 J can be as high as 42meV for SIb at 
X — 1/4, which corresponds to a Curie temperature esti- 
mated with the mean field approximation (Tcmf) of 80K 
and 44nieV for On at x = 1/9, to a Tcmf of 83K. Note 
that these ferromagnetic systems consist of no magnetic 
elements like Fe, Co, Ni, or rare-earth elements. 

The most particular feature takes place at x — 1/4 
of On when a finite MEfm is accompanied by a zero 
MEafm- That is, the initial antiparallel configurations 
in calculating Eafm were found to relax to the non-spin- 
polarized state. The zero MEafm together with a finite 
MEfm suggests the formation of polarized itinerant elec- 
trons in On, with an exchange split of 0.75eV in bands. 

The different characteristics of On and SIb at a; = 1/4 
are implied in the distinct band structures generated by 
the non-spin-polarized calculations. The band structure 
of hexagonal BN sheet has been shown previously [22I] to 
display a nearly-free-electron like characteristic around 
r point of the conduction band (CB). Since the most 
important contribution to determine the magnetic prop- 
erties in On and SIb comes from the electronic states 
near the Fermi level (Ef), we present in FIGl2]the cor- 
responding CBs for the BN sheet, and those for On and 
SIb at x — 1/4. In the BN sheet, only s-electrons of N 
atoms contribute to the free-electron like s-band (blue 
diamonds), while the p^-band (red dots) is mainly from 
the Pz -electrons of B atoms. When the O or Si substi- 
tute atoms are introduced into the system, the original 
s-band and pz-band receive contributions from both B 
and O atoms in On but from Si atom alone in SIb. Rel- 
ative to the generally inert a-band (whose energy at F 
point is taken as the energy reference zero in FIGH]), the 
effect of O and Si substitute atoms are primarily the shift 
of s- and p^-band and the change in dispersion of the Pz- 
band. In On the upward shifted s-band and the down- 
ward shifted Pz-band results in the joint of the two bands 
around the F point. The joint along with a relatively flat- 
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FIG. 2: Band structures for (a) BN sheet and (b) On (c) Sis calculated at x — 1/4. The energy zero is taken to be the 
minimum of the band a's (green diamonds). 



ter pz-hand provide an excellent condition for high DOS 
at Ep which further leads to the formation of polarized 
itinerant electrons. In SIb, the main effect is the pressed 
downward Pz-hand to below the almost unaltered s-band 
with, contrary to that in On, a more dispersive p^-band, 
i.e. a wider band width for the occupied electrons in the 
CBs. In the previous GW band structures calculations 

I 



The establishment of the polarized itinerant electrons 
in On at a; = 1/4 is further supported through the charge 
distribution for the states in CBs just below Ep around 
r point in the non- spin-polarized calculation. In On, the 
charges distribute mainly on the B atoms surrounding 
the O impurity (a rather extending distribution) which 
is further joined by the charges distributing above the 
more distant N atoms. Therefore the charges distribute 
as connected charge sheet throughout the system for the 
isosurface of density less than 0.324e[23|. This result con- 
trasts to that in the SIb system whose charges distribute 
mostly around the substitute Si atoms (a more local dis- 
tribution) and the connected isosurface occur only when 
the density is less than 0.213e. 

The x-dependence ME's are also very different for SIb 
and On- In SIb, both the MEpM and MEafm decreases 
monotonically with respect to increasing x (decreasing 
dnn)- However, in On, the two x-dependence ME's di- 
verge at values of x larger than 1/16, i.e. the decreas- 
ing and approaching zero MEafm versus an increasing 
MEpM- We should emphasize here that the behavior of 
MEpM for On differs not only from that for SIb but also 
from that for the systems carrying localized moments 
formed by other non-magnetic substitute atoms in the 
BN sheet. Examples are the 2D BN sheet with carbon 
atoms substituting for either B or N atoms. Si atoms for 
N atoms, as well as the vacancies created by removing 
either B or N atoms fToJ. The exceptional behavior of 
MEpM together with the vanishing MEafm in On pro- 



for BN sheet [2^ and h-BN[2^, the effect of many-body 
correction to the bottom of CBs is mainly in shifting the 
band energy rather than their curvatures. These results 
would support our above conclusions as the join/disjoin 
of the Pz and s-band as well as the curvature of the Pz- 
band dominate the crucial physics of magnetism in On 
and SIb. 

I 

vidcs a prominent feature for this 2D system as develop- 
ing into the polarized itinerant electron phase. 



IV. ROBUSTNESS OF THE POLARIZED 
ITINERANT ELECTRONS IN On 

In order to examine how the establishment of the po- 
larized itinerant electrons in On at a; — 1/4 depends 
on the oxygen distributions, different supercells consist- 
ing of O atoms with the same nearest-neighbor distance 
i^nn~5A but different numbers and distributions are stud- 
ied. They include a) a 6x6 supercell but with only eight 
instead of nine O substitute atoms, b) a 4x4 supercell but 
with only three instead of four O atoms c) a rectangular 
instead of a 4x2 hexagonal supercell with two O substi- 
tute atoms. All these examined systems were found to 
be still stabilized in the polarized itinerant electron phase 
with less than 4meV changes in MEfm- Regarding the 
substitute atom concentration x, the highest x of mag- 
netic On system we considered is 1/3, i.e. dnn=4.3A. At 
this concentration, the system remains stabilized in the 
polarized itinerant electron phase. However, we should 
emphasize that the role of N atoms in On is essential as 
magnetism would disappear if all N atoms are replaced 
by O atoms, i.e. the BO sheet is in fact a normal metal. 
When mapping into the 2D HEG, this corresponds to the 
paramagnetic fluid at high electron density. 

FIG [3] shows the x-dependence formation energy for 
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FIG. 3: The x-dependence formation energy for On and Sis, 
where the reference zero is chosen as the formation energy at 
X = 1/64, i.e. dnn=20A. 



On and SIb, which is defined as 
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where EoN(siB) ^^.d EBNshoct are the energies for the On 
(SiB) system (ferromagneticahy polarized ones) and the 
BN sheet calculated by using exactly the same supercell 
to eliminate numerical errors, En (Eb) and Eq (Esi) are 
the energies for the removed N (B) and the added O 
(Si) atoms respectively [l^. The formation energy for the 
system oi x — 1/64, i.e. the system with no interaction 
among the substitute atoms, is taken as the reference 
zero in the figure. Unlike the repulsive interactions in SIb 
at a; < 1/9, the O atoms introduced for substituting N 



atoms in the BN sheet exhibit attractive interaction once 
X becomes larger than 1/32. Besides, by employing a very 
large supercell (11x10) with only two substitute atoms 
(SA) created at the distance dsA-SA in the cell, it is 
shown (see the insert in FIGl3]) that for both systems the 
substitute atoms tend to aggregate to form the system of 
X = 1/4. 

Considering the possible experimental fabrication of 
this 2D system supported on a substrate, our study show 
that the system with polarized itinerant electrons can 
sustain an external applied tensile and compressive strain 
of up to 0.1, i.e. a change of 10% in lattice constant. 
Similar study for SIb, on the contrary, leads to a non- 
magnetic state when the system is under compressive 
strain. 

V. CONCLUSION 

In summary, we have explored the possible magnetic 
phases in the On and SIb systems using first-principles 
methods. The polarized itinerant electron phase was 
found in the BN3/4O1/4 system while in the SIb system 
only the paramagnetic phase at low silicon concentrations 
and ferromagnetic phase at high silicon concentrations 
were identified. The different magnetic properties in On 
and SIb were discussed through the ^-dependence mag- 
netization energy, the band structures, and the charge- 
density distributions. The robustness of the formation of 
the polarized itinerant electron phase was also discussed 
with respect to the O substitute atom distributions and 
the applied strain to the system. 
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